Symmetric correlations as seen at RHIC 
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Abstract 



We analyze the forward-backward multiplicity correlation coefficient as 
measured by STAR. We show that in the most central Au+Au collisions bins 
located symmetrically around r\ = with large separation in pseudorapidity 
O ; are more strongly correlated than bins located asymmetrically with smaller 

separation. In proton-proton collisions the opposite effect is observed. It 
suggests a qualitatively different behavior of the two-particle correlation as a 
. function of pseudorapidity sum in p+p and Au+Au collisions. 

> 
(N 
00 

1. Correlations between particles produced in different rapidity regions have 
been intensively studied since the early times of high-energy physics [1] . Particularly 
interesting are correlations between particles with large separation in rapidity. It is 
recognized that such correlations are born immediately after the collision, when the 
produced system is very small (spatial size of the order of a few femtometers) and 
! before rapid longitudinal expansion. 

One popular method to study long-range correlations is to measure the multi- 
plicity correlation coefficient, i.e., to quantify how multiplicity (number of parti- 
cles) in one rapidity window influences multiplicity in another one. This problem 
was thoroughly studied in hadron-hadron collisions at various energies [21 El E], 
[3 El El [81 [91 Uni HU H2l H31 HI]. One important lesson from these studies is that the 
forward-backward correlation coefficient decreases as a function of rapidity distance 
between bins. 

Recently the STAR Collaboration at RHIC announced the results [15] of the 
forward-backward multiplicity correlation coefficient measured in Au+Au collisions 
at yfs = 200 GeV. The measurement was performed for two narrow pseudorapidity 
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bins with the distances between them ranging from 0.2 to 1.8, covering a substan- 
tial part of the midrapidity region. For the first time very interesting features were 
observed: (i) the correlation coefficient increases significantly with centrality of the 
collision, and (ii) it remains approximately constant (except for very peripheral 
collisions) across the measured midrapidity region \r)\ < 1. These results were in- 
terpreted in the framework of the color glass condensate [IB] or the dual parton [5] 
models. 

Recently various mechanisms have been proposed to understand the data quan- 
titatively [T71 HBJ [13 ED] • However, in these calculations the sophistication of the 
STAR analysis was not fully appreciated, and the published results cannot be di- 
rectly compared with data. As emphasized by Lappi and McLerran [21] in the 
STAR analysis, the correlation coefficient is measured at a given number of parti- 
cles in an additional reference window. This procedure significantly influences the 
forward-backward correlations, and we come back to this problem later. 

In the present paper we analyze the STAR data and extend the discussion ini- 
tiated in Ref. [2T]. We describe the STAR analysis in detail and derive a general 
formula that relates the correlation coefficients measured with and without the step 
of fixing particle number in the reference window. 

The main result of this study is the observation that the two-particle pseudo- 
rapidity correlation function is qualitatively different in p+p and central Au+Au 
collisions when studied as a function of pseudorapidity sum r/i + r] 2 . In a model 
independent way we show that bins located asymmetrically around rj = with a 
small separation in pseudorapidity are significantly more weakly correlated than bins 
located symmetrically with much larger separation. It is the first time this effect is 
observed. In p+p collisions the opposite effect is observed, i.e., bins with smaller 
separation are more strongly correlated even if they are asymmetric. 

2. The multiplicity correlation coefficient for two bins X and Y is 
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d xy = {n x n Y ) - (n x ) (n Y ) ] D YY = (n Y ) - (n Y ) 2 , (2) 

where nx and n Y , respectively, are event-by-event multiplicities in X and Y. Due 
to the Cauchy-Schwarz inequality bxY varies from —1 to +1. 

The STAR Collaboration measured the multiplicity correlation coefficient be- 
tween two symmetric (with respect to n = in the center-of-mass frame) pseudora- 
pidity bins B (backward) and F (forward) of width 0.2. To reduce a trivial source 
of correlations coming from the impact parameter fluctuations^] STAR introduced 



1 Higher ub triggers a smaller impact parameter that leads to higher %. 
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the third symmetric reference bin R (see Fig. [TJ), and all averages (nB) nR , ( n B) 



'n R - 



and (nBnp) nR were measured at a given number of particles hr in this bin. Next 



they calculated the appropriate covariance and variance in the following way: 
^bf I star = Y] P(n R ) [(n B n F ) - (n B )l R ] , 



D 2 B b\star 



= Y, nn R )[{nl) nR -{n B )l 



(3) 



where P{riR) is the multiplicity distribution in the reference bin R at a given cen- 
trality class that is defined by a range of ur, i.e., rii <ur < n 2 . Equation ([3]) allows 
us to calculate the correlation coefficient as measured by STAR: 
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It is important to emphasize that if (tib), (n 2 B ), and (n B n F ) are measured without 
the step of fixing ur (namely all events are taken to directly measure D 2 BF and D 2 BB 
with hr in a given centrality range) different results are obtained^ In the following 
all observables without a label STAR denote that D BF and D 2 BB are calculated 
without fixing tir. 
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Figure 1: Configuration with maximum pseudorapidity gap between B and F. 

The STAR procedure of measuring b B p\sTAR substantially removes the impact 
parameter fluctuations, indeed. However, as shown in Ref. [21], it complicates the 
interpretation of BbfIstar since it clearly depends (in the nontrivial way) on cor- 
relations between B(F) and R. In the following we derive the relation between 
^bfIstar and multiplicity correlations bsF and bsR, = b FR that are obtained in the 
same centrality class but without the step of fixing tir. Such calculation was per- 
formed in Ref. [21J, where for simplicity the multiplicity distribution P(nB,np,nR) 
was assumed to be in a Gaussian form. Here we show that the result derived in Ref. 

2 Naively, it seems that both procedures should lead to the same result. We can always measure 
(0) nR at a given tir and calculate (O) = J2n R P{ n R) (0)n R - ^ n ^ s case j 

D 2 bf = {nsnp) - (n B ) 2 = V P (n R ) {n B n F ) - (V P{n R ) (n B ) ) , 
which is clearly different from Eq. ([3]). 
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[21] is independent on P(n B , n F , n R ) provided the average number of particles in B 
at a given is a linear function of %: 

(n B ) riR = c + CxUr. (5) 
This relation is well confirmed by STAR [22] ■ It is straightforward to show that 

c = (n B ) - (n R ) 5f£, c 1 = ^. (6) 

U RR RR 

Indeed, to obtain Eq. ([§]) both sides of Eq. (jSJ) should be multiplied first by P(n R ) 
and second by P(n R )n R and summed over n^. Using an obvious relation 

(0) nR = pT^-y £ P(n B ,n F ,n R )0, (7) 

two simple equations can be derived that allow us to calculate c and c\. 
Taking Eqs. (JHJ, (EJ), and (JJJ) into account, 



RRi 



Drf\star — D BF — c^D 
D 2 bb\star = D 2 BB -cjD RR , 



where C\ is defined in (jHJ). Consequently, b BR \sTAR is given by 

x, | _ h BF ~ b 2 BR , . 

ObF\STAR — — 72 ) l y J 

where fe^F and h BR are the appropriate correlation coefficients measured without 
fixing n R . As mentioned earlier we obtain exactly the same formula as in Ref. [21j . 
It shows that Eq. (jHJ) does not depend on P(riB,nF,nR), provided the relation (j3]) 
is satisfied. 

Here point out that the interpretation of bsFlsTAR is not straightforward. For 
example, b BF \sTAR = indicates only that b B F = b 2 BR but it does not mean that 
b B F = 0. Moreover, &rf|star can be negative even if both bsw and b B R are positive. 
We conclude that the full interpretation of &bf|star is difficult without knowing b B F 
and b BR . 

In this paper we are interested in the configuration presented in Fig. (TJ where 
the distance between B and F is a maximum one, i.e., F = [0.8 < r\ < 1], B is 
symmetric with respect to 77 = 0, and i? = [—0.5 < rj < 0.5]. In this case the 
average gap between B and R is smaller by a factor of 2 than that between B and 
.F. Assuming that the two-particle correlation function depends only on \rji — ^2! 
and is not increasing as a function of \rji — 772 1 a natural ordering b BB > b B p is 



obtained, as shown explicitly in Ref. [21] ■ Consequently 

, 1 _ b BF - b 2 BR b BR - b 2 BR _ b BR 1 
^ls™-- r3 ^<- r3157 - TT ^<2, (10) 
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since b B R < 1. In the most central collisions STAR measured 6rf|star ~ 0.58, 
which violates this boundjf] Thus we arrive at an interesting conclusion that in the 
midrapidity region in the most central Au+Au collisions the following inequality 
holds: 

b B R < b BF - (11) 

It was checked by STAR that narrowing the reference bin R from \rj\ < 0.5 to 
\r]\ < 0.1 (so that all windows have the same widths) slightly increases the corre- 
lation coefficient b B F\sTAR- Also an alternative method of centrality determination 
was carried out using the STAR zero-degree calorimeter (measurement of forward 
neutrons) for the — 10% centrality, and b B F\sTAR is very close to |. In this case the 
same formula (j3J) applies; however, there are no explicate cuts on n R . We conclude 
that the width of R and the centrality cut on n R is not a factor in the result (llip . 

3. It is interesting to estimate the numerical values of the correlation coefficients 
bsF and b BR . As mentioned earlier we are mostly interested in the configuration 
where the distance between B and F is a maximum one (At/ = 1.8 in the STAR 
notation) and R is defined by |t/| < 0.5. 

As seen from Eq. ((§]) evaluation of b B p = D 2 BF /D 2 BB is straightforward. The 
covariance D bf \star and variance D bb \star are published in [15] (only for — 10% 
centrality bin). From Ref. [22j one sees that (n B ) nR is a linear function of n R with a 
coefficient C\ ~ 0.2. To calculate D 2 RR = (n 2 R ) — (n R ) 2 we use the uncorrected (raw) 
multiplicity distribution P(n r R w ) as published in Ref. [23], and take the efficiency 
correction to be n R /n T R w = 1.22 j22j [23]. Performing a straightforward calculation 
we obtain^ D 2 RR rs 4320, which allows us to calculate b B F- Taking Eq. ([9]), b B F, and 
measured b BB \sTAR into account we obtain 

b BR « 0.58, b BF « 0.72. (12) 

As seen from ( 1T2"|) in the most central Au+Au collisions b BR is significantly smaller 
than bsF- Let us note here that the average distance between B and R (one unit of 
rf) is smaller by a factor of two than that between B and F. 

It is also interesting to see how b BF depends on the distance At/ between bins B 
and F. Taking Eq. (JSJ) into account and repeating calculation^] presented above we 

3 The STAR result has an uncertainty ±0.06. Even if one assumes that the measured bsFlsTAR 
is slightly below 0.5, it is still difficult to understand with an assumption bsR > bsF, since it 
requires bsR ~ bsF ~ 1. 

4 We take P(n™ w ) oc exp(~^) for 431 < ng w < 560 and P(n™ w ) oc exp(- {n \^ 1)2 ) for 
n™ w > 561, which gives D 2 RR \ raw = 2904. Consequently, D 2 RR = (l.22 2 )D 2 RR \ Taw . 

5 For small Aij the reference window R is composed of two windows 0.5 < \r]\ < 1 and we assume 
that c 2 D RR is approximately the same as with R defined by < 0.5. 
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found that bsF in central Au+Au collisions is approximately constant as a function 
of Arj, which is consistent with the dependence of &bf|star on Arj. 

Finally, let us notice that STAR also measured bs f \ star in P+P collisions; how- 
ever, in this case the exact value of C\ is not known. We checked that for a very 
broad range of C\ we always obtain a standard ordering b BB > &bf@ 

4. Several comments are warranted: 

(i) To calculate the correlation coefficients bsF and bsR the experimental values 
of D 2 bf \star and D bb \star are required as an input. Unfortunately they are pro- 
vided only for the most central collisions. It would be interesting to measure the 
centrality dependence of the effect reported in this paper. It is expected that in 
peripheral collisions the standard relation b B R > bsF should be recovered. If so, it 
would indicate a qualitatively different behavior of central and peripheral Au+Au 
collisions. 

(ii) It is worth mentioning that HIJING [H] and the Parton String Model (PSM) 
[25J fail to describe the Au+Au data for the forward-backward multiplicity correla- 
tion coefficient. However, they are consistent with the p+p data. In the most central 
Au+Au collisions, and for the configuration presented in Fig. [U both models predict 
bBFlsTAR < \i which is consistent with the relation bsR > b BF ^\ 

(iii) It is not straightforward to propose a realistic mechanism that more strongly 
correlates bins B and F than bins B and R. One possible mechanism is the formation 
of certain clusters strongly peaked at rj = that decay symmetrically into two 
particles. This mechanism obviously correlates bins B and F and introduces no (or 
much weaker) correlations between bins B and R. To go beyond speculations more 
detailed measurement of the forward-backward correlations between symmetric and 
asymmetric bins is warranted. 

5. In summary, we analyzed the STAR data on the forward-backward multi- 
plicity correlation coefficient &bf|star in the most central Au+Au collisions. This 
measurement was performed with the intermediate step of fixing the number of 
particles in the third reference window R, see Fig. (TJ and we emphasized the im- 
portance of this step. We derived the general formula that relates &bf|star and the 
correlation coefficients bsF and bsR measured in B — F and B — R without fixing 
the number of particles in R. 

The most important result is the observation that for the configuration presented 
in Fig. [T| in the most central Au+Au collisions, the correlation coefficient bsR is 

6 We assume P(n#) to be given by a negative binomial distribution with standard parameters 
(nn) = 2.3 and k = 2. Taking, e.g., c\ = 0.1 we obtain bsR ~ 0.28 and bsF ~ 0.13. 

7 In particular b BF \ ST AB. ~ 0.1 in HIJING and b BF \ STAR sa 0.4 in PSM, see Ref. [TS]. 
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significantly smaller than bsF- This is exactly opposite of what is expected and 
measured in p+p collisions (the distance between B and R is smaller by a factor of 
2 than that between B and F). Moreover, we found that in central Au+Au colli- 
sions, bsF is approximately constant as a function of the pseudorapidity separation 
between symmetrically located bins B and F. To understand these results it is 
necessary to assume that in central Au+Au collisions the two-particle correlation 
function strongly decreases as a function of 1 771 + 772 1 - It indicates the presence of 
a specific mechanism of correlation that strongly correlates bins located symmetri- 
cally around 77 = for which \r]i + 772 1 ~ 0, but is less effective for asymmetric bins 
\Vi + V2\ > 00 

In this paper we solely concentrated on an analysis of the experimental results 
and at the moment we see no compelling explanation of this effect. It would be 
interesting to directly measure at RHIC and LHC the multiplicity correlation coef- 
ficient for symmetric and asymmetric bins to confirm conclusions presented in this 
paper. 
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